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Abstract 

Very light right-handed (RH) sneutrinos in the Next-to-Minimal Supersym- 
metric Standard Model can be viable candidates for cold dark matter. We in- 
vestigate the prospects for their direct detection, addressing their compatibility 
with the recent signal observed by the CoGeNT detector, and study the impli- 
cations for Higgs phenomenology. We find that in order to reproduce the correct 
relic abundance very light RH sneutrinos can annihilate into either a fermion- 
antifermion pair, very light pseudoscalar Higgses or RH neutrinos. If the main 
annihilation channel is into fermions, we point out that RH sneutrinos could 
naturally account for the CoGeNT signal. Furthermore, the lightest Higgs has a 
very large invisible decay width, and in some cases the second-lightest Higgs too. 
On the other hand, if the RH sneutrino annihilates mostly into pseudoscalars or 
RH neutrinos the predictions for direct detection are below the current experi- 
mental sensitivities and satisfy the constraints set by CDMS and XENON. We 
also calculate the gamma ray flux from RH sneutrino annihilation in the Galactic 
centre, including as an interesting new possibility RH neutrinos in the flnal state. 
These are produced through a resonance with the Higgs and the resulting flux 
can exhibit a significant Breit-Wigner enhancement. 



*MultiDark Fellow 



1 Introduction 



Very light weakly-interacting massive particles (WIMPs) are currently receiving much 
attention as a potential solution to the dark matter problem. This has been motivated 
by some recent experimental results in direct detection experiments that might favour 
them over more conventional scenarios with heavier candidates. 

In their search for the elastic scattering of dark matter particles, the CoGeNT 
collaboration observed an irreducible excess in their data [1] that, if interpreted in terms 
of WIMPs, would correspond to a very light particle, with mass in the range 7—12 GeV, 
and a large elastic scattering cross section, of order 10~^ pb. Furthermore, hints of an 
annual modulation in the CoGeNT experiment have also been observed [2], which 
are also consistent with their previous results and narrow the range of WIMP masses 
down to 7 — 9 GeV. Although, such a particle has similar properties to the candidate 
suggested to account for the annual modulation signal reported by the DAMA/LIBRA 
collaboration [3, 4], a joint explanation of both experimental results is not possible 
unless extreme assumptions are made for the different uncertainties [5, 6, 7, 8] such as 
the inclusion of large quenching factors or channeling effects. 

These observations are, however, challenged by the negative results obtained in 
searches by other experimental collaborations. Most notably, CDMS [9], XENONIO 
[10], XENONIOO [11] and recently, SIMPLE [12] have set upper bounds on the spin- 
independent part of the WIMP-proton cross section that are in tension with the regions 
of the parameter space compatible with DAMA/LIBRA and CoGeNT. The compat- 
ibility between these experimental results was tested in Ref. [13] using Bayesian sta- 
tistical methods. It was claimed that when uncertainties in the scintillation efficiency 
of XENONIOO are taken into account, the resulting exclusion limit is not sufficient 
to rule out the CoGeNT region (see also Ref. [14]). Also, from a theoretical point of 
view, an unconventional dark matter candidate, coupling differently to protons and 
neutrons, could account for the CoGeNT signal while having escaped detection in 
XENON [15, 16]. Regarding the comparison of CoGeNT with CDMS data, there is no 
consensus between both collaborations on to which extent their spectra for low-energy 
events observed in both experiments disagree, see in this respect Refs. [17] and [18]. Al- 
though channeling effects in the CoGeNT crystals could help reconciling both results, 
it is not clear if this effect can be large enough [19]. 

Various theoretical constructions with very light WIMP dark matter have been 
proposed in the literature. For example, in the case of supersymmetric models, it 
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was realised that very light neutralinos were viable in the Minimal Supersymmetric 
Standard Model (MSSM) both in an effective low-energy description [20, 21, 22, 23, 24] 
as well as with parameters described at the Grand Unification scale (see e.g., [25]), 
provided the mass terms for both scalars and gauginos were no n- universal. These 
neutralinos can lie within the DAMA/LIBRA region [21, 22, 26, 27] when the correct 
relic density is required. However, this possibility relies on the use of light scalar and 
pseudoscalar Higgses in order to enhance the neutralino annihilation cross section, and 
a choice of parameters that leads to sizable contributions to low-energy observables. 
In particular, the branching ratios of some rare decays (mainly 6 — > 57 and Bs — )■ 
H~^fi^) impose very stringent constraints that make these solutions very fine-tuned [28]. 
Furthermore, large values of the ratio of the Higgs vacuum expectation values (tan /3) 
are becoming increasingly constrained by the recent results from the LHC [29, 30]. 

This problem can be in principle alleviated when extensions of the MSSM are 
considered. In particular, in the Next-to-MSSM, where an extra singlet field is included 
in order to provide a fi parameter of order of the Electroweak scale, the contribution to 
low-energy observables can be reduced. Very light neutralinos are also possible within 
the NMSSM [31]. An increase in their annihilation cross section can be obtained in 
the presence of either a very light scalar or pseudoscalar Higgs. These can be viable 
if their singlet composition is large enough. The predictions for the neutralino-nucleus 
scattering cross section, which are very sensitive to changes in the Higgs sector [32, 33], 
vary significantly and span several orders of magnitude in the small mass region [34]. 
Recent analyses of this scenario have been made [35, 36, 37, 38, 39, 40] in which the role 
of experimental constraints and the naturalness of the parameters have been thoroughly 
addressed in the light of the CoGeNT results. 

The neutralino is not the only viable WIMP in supersymmetric theories. Another 
interesting possibility is the sneutrino [41]. The sneutrino is excluded as dark matter 
in the MSSM, where it is a left-handed (LH) field, due to its large annihilation cross 
section and excessive scattering cross section off nuclei [42]. Nevertheless, it can be 
viable in extended models where its coupling to the Z boson is reduced. This can be 
achieved by introducing a right-handed (RH) sneutrino superfield (which also entails 
the inclusion of a RH neutrino and thus the possibility of having a see-saw mechanism 
that accounts for the smallness of neutrino masses). The sneutrino can thus be a mixed 
LH-RH field [43, 44, 45, 46] and reproduce the correct relic density, even if it is very 
light [47]. However this requires a particular supersymmetry-breaking scheme with very 
large trilinear terms which is not available with the standard supergravity mediated 
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paradigm. Another possibility is having a pure RH sneutrino, which is generically 
non-thermal due to its extremely small Yukawa coupling [48, 49, 50, 51] unless it is 
somehow coupled to the observable sector, for example via an extension of the gauge 
[52, 53, 54, 55] or Higgs [56, 57, 58, 59] sectors. 

In this paper we will work with an extension of the MSSM that was presented in 
Ref. [57] and in which two new singlet superfields were included, as in Refs. [60, 56]. An 
extra singlet superfield 5* addresses the /i problem in the same way as in the NMSSM 
and provides extra Higgs and neutralino states, while an extra singlet superfield N 
accounts for RH neutrino and sneutrino states. The phenomenology of this construction 
was studied in [59], where the possibility for very light RH sneutrinos as candidates for 
cold dark matter was already pointed out. 

We will further investigate here the viability of very light RH sneutrinos in the light 
of the current experimental situation. In particular, in Section 2 we will review the 
conditions under which the correct sneutrino relic density is achieved. We find three 
possible scenarios, where RH sneutrinos annihilate mainly in a fermion-antifermion 
pair, very light pseudoscalar Higgses or RH neutrinos. We explore these annihilation 
channels and investigate the available parameter space in each case. Then, in Section 3 
we compute the spin-independent part of the RH sneutrino-nucleon elastic scatter- 
ing cross section. In order to contemplate all the possibilities, we will consider cases 
where the RH sneutrinos can have an elastic cross section off nuclei which is large 
enough to provide a WIMP interpretation of the CoGeNT signal but we will also study 
other scenarios in which the predictions are below the current exclusion limits set by 
CDMS and XENON. In an attempt to find discriminating features of this scenario 
that distinguishes it from the case of very light neutralinos we will study some possible 
collider implications in Section 4, computing the invisible decay width of the lightest 
and second-lightest CP-even Higgs bosons. After that, in Section 5 we compute the 
gamma ray flux from RH sneutrino annihilation in the Galactic Centre (GC) for each 
of the scenarios, and compare it to the current data by the Fermi satellite. As a novel 
feature, the channel with RH neutrinos in the final state is included in the analysis. 
Finally, our conclusions will be presented in Section 6. 
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2 Very light RH sneutrinos in the NMSSM 



The details of the model with RH sneutrinos in the NMSSM were already introduced 
in Ref. [57] and the phenomenology of this supersymmetric WIMP was studied in more 
detail by some of us in Ref. [59]. In particular, we showed that the correct relic density 
could be obtained for a wide range of sneutrino masses and in extensive regions of the 
NMSSM parameter space. 

The superpotential of this construction is expressed as 

W = Wt^MSSM + XnSNN + yNL- H2N, 
M^NMSSM = Y^H2-Qu + YaH^-Qd + Y,H,-Le-XSH,-H2 + ^KS\ (2.1) 

where flavour indices are omitted and the dot denotes the SU (2) l antisymmetric prod- 
uct. As in the NMSSM, a global Z3 symmetry is imposed for each superfield, so that 
there are no supersymmetric mass terms in the superpotential. After radiative Elec- 
troweak symmetry-breaking the Higgs fields get non-vanishing vacuum expectation 
values (VEV). The VEV of the singlet Higgs, Vg, induces an effective /i parameter, 
H = Xvg, and a Majorana mass for RH neutrinos, = 2Xi^Vs, both of order of the 
Electroweak scale. The phenomenology of this model is largely dependent on the SNN 
coupling and thus on the new parameter Ajy. 

The Lagrangian, with the corresponding soft-supersymmetry breaking terms reads 

r 2\r\\2 I 2|~|2 I 2|J|2 I 2|rl2 , 2|~|2 

-Acalarmass = "^qIQI + +"^^1^1 +"^e|e| 

+m^^ \H^f + m\^\E2f + m\\S\^ + m\\N\\ (2.2) 
-£A-terms = [a^YuH^Qu + AdYdHiQd + AeYMiLc + H.c.) 

+ (^-XAxSH^H^ + ^t^A^S"" + H.c.^ 

-I- (A^v^A^^iV^ + y^Ay^LH^N + H.c.) , (2.3) 
and contains a soft mass term for the RH sneutrino, m^y, and new trilinear soft terms 

■^Xn ^^'^ ^VN- 

The smallness of the neutrino Yukawa coupling implies that the sneutrino mass 
eigenstates have a negligible mixing and therefore can be identified with the LH and 
RH components. The lightest sneutrino is then a pure RH field and its mass can be 
expressed in terms of the NMSSM parameters as follows 

= + |2^7v^sP + \yNV2\^ + 2AAr {Axj^v^ + (kt;^ - XviV2y) . (2.4) 
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Figure 1: Trajectories in the (m^,AAr) plane with fixed RH sneutrino mass, given various 
values of Axj^. For each choice of Axj^ the dashed line represents the trajectory along 
which = 50 GeV and the solid one corresponds to m^^ = 0. We have used tan/3 = 5, 
A = 0.3, K = 0.2, and /i = 200 GeV. 

The flexibility of this construction stems from the fact that the new free parameters 
{Xn, ^Ajv) can be chosen in such a way that they provide a wide range of RH 
sneutrino masses and couplings while on the other hand not affecting the rest of the 
NMSSM spectrum. This is illustrated in Fig. 1, where we show the trajectories for a 
fixed sneutrino mass in the (m^, Ajv) plane for various choices of the trilinear parameter. 
We observe that very small RH sneutrino masses can be obtained for any choice of the 
soft mass parameter (for example, increasing \Axj^\ larger values of Aa? are possible). 
This is potentially interesting, since it suggests that no large non-universalities in 
the soft parameters are needed in order to have a very light RH sneutrino (contrary 
to what happens with very light neutralinos, where the bino mass parameter has to 
be significantly reduced). For concreteness, the soft masses for sleptons are fixed to 
nT'L,R = 250 GeV, squark masses are set to mQ^u,D = 1 TeV or 1.3 TeV, depending on 
the specific example, although our conclusions regarding dark matter are not sensitive 
to these choices. We also adopt the GUT relation for gaugino masses at low energy. 
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In order for very light RH sneutrinos to be viable dark matter candidates they 
have to reproduce the correct value for their thermal relic abundance, Qf^h"^. In [59] 
some of us showed that there were three distinct scenarios in which sneutrinos with 
masses below 10 GeV could be in agreement with the constraint set (at the 2 a level) 
by the WMAP satellite 0.1008 < Vlffh"^ < 0.1232 [61]. The three cases correspond to 
RH sneutrinos annihilating preferentially either in fermions (A^A^ — )■ //, mainly into 
bb), or in a pair of very light pseudoscalar Higgses (A^A^ — )■ A^A^), or in RH neutrinos 
{NN NN). 

In our calculations the analysis of the NMSSM phenomenology has been performed 
with the NMHDECAY 2.3.7 code [62, 63] , which minimises the scalar potential, dismissing 
the pressence of tachyons and/or false minima, and computes the Higgs boson masses 
including 1- and 2-loop radiative corrections, as well as the rest of the supersymmetric 
masses. It also implements the different collider constraints that apply to the Higgs 
sector. Based on this code, we have built a set of routines which numerically calculate 
the RH sneutrino spectrum and relic density, which is then compared to the obser- 
vational bound extracted from the WMAP results. For details of our calculation we 
refer the reader to Ref. [59], where the amplitudes for RH sneutrino annihilation are 
explicitly computed for each possible channel. Experimental constraints on low energy 
observables are also taken into account. In particular we incorporate the recent bounds 
on rare decays {b — > S7 and Bs f^~^f^ ) and the muon anomalous magnetic moment. 

In particular, we impose the experimental bound on the branching ratio of the rare 
b s-f decay, 2.85 x 10^"^ < BR{b 57) < 4.25 x 10"^ at 2a level, obtained from 
the experimental world average reported by the Heavy Flavour Averaging Group [64], 
and the theoretical calculation in the Standard Model [65, 66], with errors combined 
in quadrature. We also take into account the recent upper constraint on the Bs 

branching ratio obtained by the CMS collaboration, BR(i?5 — )■ yU+yU^) < 1.9 x 
10~® at 95% c.l. [67], which improves the previous results from the LHCb and CDF 
collaborations. 

Concerning the muon anomalous magnetic moment, a constraint on the supersym- 
metric contribution to this observable, a^^^^, can be extracted by comparing the exper- 
imental result [68], with the theoretical evaluations of the Standard Model contributions 
[69]. When e~^e~ data are used the experimental excess in = {g^ — 2)/2 would con- 
strain a possible supersymmetric contribution to be 3.3 x 10~^*^ < a^^^ < 47.9 x 10^^" 
at the 2a level, where theoretical and experimental errors have been combined in 
quadrature. However, when tau data are used a smaller discrepancy with the exper- 
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Table 1: Sets of inputs corresponding to the examples used in the analysis. The resulting 
masses of the CP-even and CP-odd Higgses are indicated, together with the composition 
of the CP-even Higgses. All the masses are in GeV. The corresponding values for low 
energy observables are also indicated. In scans nnl) and nn2) the value of A is determined 
by requiring a specific RH neutrino mass as a function of the rest of the parameters, 
Mat = 8 GeV and = 15 GeV in examples nnl) and nn2), respectively. See the text for 
more details. 
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imental measurement is found [70]. Due to this reason, in our analysis we will not 
impose this constraint, but only indicate the resulting prediction in our benchmark 
examples. 

Finally, the presence of an new unrealistic vacuum in the NMSSM field space which 
is not included in the NMHDECAY code has been pointed out in Ref. [71]. It was found 
that the direction where {Hi) = {H2) 7^ and (S) 7^ 0, with vanishing D-terms, is 
equivalent to the so-called MSSM unbounded from below (UFB)-l direction [72], but 
lifted up by the singlet Higgs S in the NMSSM. This means that the UFB-1 in the 
MSSM becomes just another minimum in the NMSSM. Imposing that the realistic 
minimum is deeper than this new one leads to constraints that can be very stringent 
for large values of the Ax and parameters. In our analysis we have implemented 
these constraints too. 

Let us now explain in more detail how the correct annihilation cross section can 
be obtained for very light RH sneutrinos. In Table 1 we detail the inputs of a set of 
benchmark scenarios that will be used throughout the text. We also include information 
about the resulting Higgs spectrum. 

2.1 NN ff 

This case is very simple to analyse. There is only one Feynman diagram that con- 
tributes, namely the exchange of a CP-even Higgs, along the s-channel depicted 
in Fig. 2 (the s-channel annihilation mediated by the Z boson vanishes for a pure RH 
sneutrino). For light sneutrinos the main annihilation product is a bb pair. Annihila- 
tion into cc can also be significant, as we discuss later, depending on the composition 
of the lightest Higgs. 

Under these circumstances it is easy to derive an analytical approximate expression 
for the thermally averaged annihilation cross section by means of a partial wave expan- 
sion (avMoi) ~ a + bx, where x = T/m is proportional to the WIMP velocity-square. 
This approximation holds when one is far enough from resonances and thresholds for 
new final states^. 

For this specific diagram the integral of the matrix element describing each annihi- 
lation process, A^A^ — )■ X1X2, which we define in terms of the scattering angle in the 

^Although for very light sneutrinos resonances can generally be avoided (as they would require very 
light CP-even Higgses), the threshold for annihilation into bb happens around 4 GeV and therefore 
we should expect deviations from this approximation for sneutrinos lighter than 6 GeV. 
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Figure 2: Diagram contributing to the annihilation of RH sneutrinos into //. 

CM frame, 6cm ^ as 

wx.xM = \ COS OcM \AiNN ^ XiX2)p, (2.5) 



can be written as 



grrih V '^HONN'^HONN^Hf^HO 

2Mwcosf3 J \j 



— ^ {2s - 8ml) , (2.6) 



Wr 



2^ X-. - ^^c) , (2-7) 



2Mw sin /3 / A, 



where C^ofjf; is the RH sneutrino couphng to the Higgs and Ajj is the square 
denominator of the Higgs propagator, both defined in Appendix A of Ref. [59]. The 
following convention is used to express the composition of the CP-even Higgs mass 
eigenstates, i/f = 5^0 Hd + 5*^0 + 5^0 S. 

i i i 

In terms of the quantities 1^X1X2 the annihilation cross section can be calculated 
numerically as detailed in Appendix B of Ref. [59]. Moreover, the coefficients of the 
partial wave expansion of the thermally averaged annihilation cross section can also be 
calculated analytically [73]. Following the prescription of Ref. [74], the expressions for 
the velocity-independent contribution to the annihilation cross section into a pair of hh 
or cc then read 



47r V2MvKCOs/3y 

2 [im?- — m^"!^/^ 



47r V2MvKsin/3y ' ^ 



where we have defined 



V ; ^ = V- 2.9 
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In most cases the lightest Higgs contribution will dominate in the expressions above 
(especially if its mass is small), so that we can define 

cl q2 

V^V^^ Y I ; U^U,^ ^ Y . (2.10) 

The contribution from the h parameter in the {(JVmoi) expansion is suppressed (since 
~ 1/20) and will be neglected in the following discussion. 

Notice that 



2 _2\ 3/2 



for RH sneutrinos with a mass ^7 — 9 GeV. Thus, annihilation into cc can only 
be dominant for small values of tan /3 and depending on the properties of the Higgs 
sector. Since it is the lightest Higgs the one that contributes the most to both V and 
U, the above condition could happen if its Hu component is much larger than its Hd 
component. Although this is not a very common situation in our scans, we will keep 
this possibility open and explicitly consider both cases, where annihilation into either 
bb or cc dominates^. 

For moderately heavy WIMPs the sneutrino relic density is usually approximated 

as 

' L07xlO»GeV-' 

where Mp = 1.22 x 10^^ GeV is the Planck mass and g*{xf) is the number of relativistic 
degrees of freedom at the decoupling temperature. However, very light dark matter 
with mass smaller or of order 10 GeV would have decoupled when the temperature 
of the Universe was around 400 MeV, precisely when quarks confined into hadrons. 
This hadronization implies that the number of relativistic degrees of freedom drops 
dramatically around the decoupling temperature of very light WIMPs, producing an 
enhancement of their relic abundance [75]. Furthermore, this enhancement is sensitive 
to the phase transition model that is considered (in particular, it depends on the 
deconfinement temperature, Tc). In our analysis we have taken Tc = 400 MeV. 



^It should be noted that the contribution from annihilation into tt can be larger than that cor- 
responding to cc. However, it has the same dependence on the lightest Higgs composition than the 
contribution from b quarks (i.e., Orr is proportional to D^) and because of the different Yukawa cou- 
plings Urf <C ttfji is always satisfied. Thus, we only have two possible regimes, where either bb or cc is 
the leading contribution. 
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For a quick analytical estimate let us momentarily assume that the QCD transition 
is close to first order and that the relativistic number of degrees of freedom suddenly 
varies from ^ 9 to 3.7, considering also the uncertainty in Tc. Notice that this 
implies essentially an increase of a factor 2.5 in the approximate evaluation of Qf^h"^ of 
Eq. (2.12). In order to reproduce the WMAP result, this results in a condition on the 
annihilation cross section that can be written as 

((Ti;) ^ a;,fe ^ 0.77- 1.9 pb if I?>-^W, (2.13) 

0.4 

(at;) ^ a^g ~ 0.77- 1.9pb if P< -U , (2.14) 

tan p 

where the lower (upper) value applies to a RH sneutrino which is heavier (lighter) than 
~ 8 GeV and therefore decouples above (below) Tc. These values are generic for 
a dark matter particle in which s-wave annihilation dominates [76] . 

So far we have used an analytical approach with several approximations so that 
the correlation between the sneutrino annihilation cross section and its scattering cross 
section off nuclei (which we calculate in the next section) is manifest. However, we stress 
that in our results we use the full calculation of VL^jh"^ following the same numerical 
method that we detail in Ref. [59]. Also, the QCD transition is taken into account 
through a parametrization of g*{xf) according to [77] with Tc = 400 MeV. 

The question is then whether or not it is possible to find a sufficiently large an- 
nihilation cross section in this model and which is the choice of input parameters to 
achieve this. The quantities a^,i and are very sensitive to the structure of the Higgs 
sector and the new couplings C^i^o^^, and our model provides much flexibility in this 
sense. For example, both a^j and Ccg increase as the mass of the lightest Higgs (which 
is the leading term) decreases, but they also depend on the composition of this lightest 
Higgs. In the NMSSM it is possible to have a very light Higgs without violating the 
current experimental limits as long as its singlet composition is large enough. 

Thus one possible scenario in which a sufficiently large annihilation cross section 
can be achieved involves a light singlet-like Higgs. The if„ or Hd components, though 
small, would determine whether predominant annihilation occurs via hb or cc, according 
to the expressions above. This kind of scenarios is interesting since the Higgs structure 
is completely different to that of the MSSM for very light neutralinos (in which Higgses 
are in the so-called intense coupling regime for which their masses are very similar and 
of order 100 GeV). On the other hand, as already commented in the Introduction, a 
light singlet-like Higgs is also one of the possible solutions for the very light neutralinos 
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in the NMSSM. 

In Table 1 we display a specific example with these properties, namely case ffl). The 
lightest Higgs, with a mass of 62.4 GeV is mostly singlet, whereas the second lightest 
Higgs is Standard Model (SM)-like with a predominant Hu composition (S'?ro = 0.98) 
and a mass of 119.4 GeV. The sneutrino relic abundance is represented on the left-hand 
side of Fig. 3 as a function of the sneutrino mass for A at = 0.25 and A\j^ = —500 GeV. 
The solid line corresponds to the estimation using the partial wave expansion of (crv). 
The dashed line corresponds to the approximation {av) ^ aj,], and the dotted line 
represents the contribution coming from only the lightest Higgs to a^,^. In extracting 
these approximations we have assumed a sharp QCD transition with Tc = 400 MeV. 
As we can see, this is one of the examples for which the contribution of the lightest 
Higgs provides a good approximation to the total result. The predictions from our 
numerical code are represented by circles. The rapid decrease of the relic density 
towards larger masses is due to the resonance with the lightest Higgs which takes place 
for mff_^ ^ mfjo/2 ^31 GeV. 

Scanning in the Xn and Ax^ parameters we can vary the Cj^of^ff coupling in such 
a way that the correct relic density is obtained for a range of RH sneutrino masses. In 
order to illustrate this, on the right-hand side of Fig. 3 we represent the distribution 
of points in the {Xn,Axj^) plane which reproduce the correct relic abundance for RH 
sneutrinos in the mass range =5 — 25 GeV. 

Another possible scenario consists of having a lightest Higgs with a mass of order 
114-120 GeV and which is mostly (as in the MSSM) and then increasing C^joffff by 
taking larger values of the parameter Aat (as we showed in Fig. 1, a larger \Ax^\ may 
then be needed in order to keep the small values of the RH sneutrino mass). 

These features are clearly seen in Fig. 4, where we consider the example ff2) in 
Table 1, which falls into this category. The relic density is plotted as a function of the 
sneutrino mass on the left-hand side of the figure. In this particular case, the masses 
of the lightest Higgs (116 GeV) and second lightest Higgs (159 GeV) are relatively 
close to each other, and both have to be taken into account when computing the relic 
density. Thus, whereas the approximation of considering only the term a^^ (dashed 
line) remains valid, the contribution from including only the lightest Higgs deviates 
significantly from the correct result (the line is not represented in the plot, as it leads 
to a relic density which is one order of magnitude above the exact value). Notice also, 
that although in this scenario the component of the lightest Higgs is not small, it 
is not sufficiently larger than the Hd component for the annihilation into cc to become 
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Figure 3: Results for case ffl) of Table 1. Left) Sneutrino relic density as a function of 
the sneutrino mass for Ajv = 0.25 and A\j^ = —500 GeV. The solid line corresponds to a 
calculation using the partial wave expansion of {crv), the dashed line uses the approximation 
(av) ~ afji and the dotted line represents the contribution coming from only the lightest 
Higgs to a^g. The circles correspond to the results using a numerical evaluation of the 
relic abundance with a parametrization of g* in which Tc = 400 MeV. Right) Values of the 
parameters and Axj^ for which the WMAP relic density is reproduced for sneutrinos in 
the mass range m^y^ =5 — 25 GeV. 

comparable to the bb channel. This would only be possible for smaller values of tan /3 
and a much more careful choice of the parameters entering the Higgs mass matrix. 

Regarding the allowed regions in the {Xn, ^Xn) plane, these are represented on the 
right-hand side of the figure for sneutrino masses in the range from 5 to 25 GeV and, 
as mentioned above, feature larger values of the A ^ parameter in order to compensate 
for the larger Higgs mass. In general it is still possible to reproduce the correct relic 
abundance for very light sneutrinos for this range of sneutrino masses and a wide area 
of the parameter space. 

Summarising, the conditions under which very light RH sneutrinos can have the 
correct relic abundance when they annihilate into a fermion-antifermion pair are rel- 
atively easy to fulfil in our model. Of course, this is due to the flexibility of the new 
input parameters Xn, nij^, and Axj^, which modify the RH sneutrino mass and cou- 
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Figure 4: The same as in Fig. 3 but for case ff2) of Table 1 with = 0.5 and Axj^ = 
-500 GeV. 

plings without altering any other feature of the NMSSM spectrum. Thus significant 
constraint appears in the new parameters. The low-energy constraints can be fixed 
only with an adequate choice of the NMSSM inputs and the RH sneutrino sector freely 
modified after that. Something that should be emphasized is that since A^r can be cho- 
sen to be rather large, the Higgs coupling to b quarks needs not be too large in order to 
reproduce the correct relic abundance, in other words, the value of tan /3 can be kept 
small. This is an important difference with respect to the case of the neutralino. 

2.2 NN 

The second possibility for very light RH sneutrinos to have the correct relic abundance 
is that they annihilate predominantly into a pair of very light pseudoscalar Higgses. 
This channel is possible provided the pseudoscalar is mostly singlet, thereby evading 
present experimental constraints and requires some tuning of the NMSSM parameters. 

In Ref. [31] it was already shown that a very light pseudoscalar Higgs could be viable 
in the NMSSM and that this made it possible for the lightest neutralino to have the 
correct relic density in this extension of the MSSM. In particular, the neutralino annihi- 
lation cross section is enhanced through a resonance with the pseudoscalar Higgs when 
2m^o ^ m^o or because the annihilation channel into A^A^ becomes kinematically 
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Figure 5: Diagrams contributing to the annihilation of RH sneutrinos into 
allowed. 

Obtaining very light pseudoscalars in the NMSSM requires a careful tuning of some 
of the parameters so that either the f/(l)i?, or U{V)pq symmetries of the model are 
quasi-restored and the CP-odd Higgs corresponds to the pseudo-Goldstone boson of 
the broken symmetry [31]. For example, this can be achieved by taking /t — or the 
trilinear terms Aa, 0. Very light pseudoscalars are very constrained by collider 

searches and its composition has to be mostly singlet-like in order to avoid these. In 
particular, they can lead to observable signals in radiative T decays. The most recent 
data from the CLEO collaboration [78] set stringent upper bounds whose effect in the 
NMSSM parameter space was investigated in Refs. [79, 80] for the case m^o < 9.2 GeV 
and later analysed in more detail for pseudoscalar masses between 9.2 and 10.5 GeV in 
[81]. The constraints coming from the analysis of Ref. [79] are incorporated in the code 
NMHDECAY 2.3.7, which we use in our analysis. The bounds from Ref. [81] are easy to 
implement in our analysis. In all cases, these are avoided if the lightest pseudoscalar 
is singlet-like. 

All the above can be applied to our model, since the inclusion of RH sneutrinos has 
no influence on any aspect of the NMSSM spectrum (or the NMSSM vacuum [59]). 
Notice however that the RH sneutrino has no coupling with the CP-odd Higgs and 
therefore there is no s-channel annihilation mediated by this particle (and therefore 
no resonant enhancement as in the case of the neutralino). It is possible, nevertheless, 
that a RH sneutrino particle annihilates preferentially into a pair A'jA'j', and this is the 
case we study here. 

The Feynman diagrams which are involved in this process are represented in Fig. 5 
and consist of a quartic coupling and s-channel mediated by CP-even Higgses. The 
explicit expression for ty^Oyio can be found in Ref. [59]. It can easily be seen that in 
the partial wave expansion of the annihilation cross section the a^io^o term is non- 
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Figure 6: The same as in the right-hand side of Fig. 3 but for cases aal) and aa2) of Table 1 
on the left and right, respectively. 

vanishing and is generally a good approximation. The relevant couplings are now 
the quartic coupling C^o^oj^j;^ and the CP-even coupling to a pair of CP-odd Higgses 
Caoa°ho- For singlet-like pseudoscalars, C'^oyio^Tv ~ ^N^t^ — 2Xn), whereas C^Oj^oho is 
much more sensitive to the CP-even Higgs composition but is independent on the Xm 
parameter. 

The contribution from the s-channel is generally sizable since it involves the VEVs 
of the scalar Higgses. For example, in the case of a pure singlet pseudoscalar and a 
lightest Higgs which is mostly one has C^o^Ofjo ~ i^Xnvi + X'^V2). This channel, if 
open, is easily more important than the NN — )■ // channel discussed in the previous 
section. Furthermore, since this coupling is more effective, the value of Atv (which now 
only appears through the sneutrino-sneutrino-Higgs vertex) that is necessary in order 
to reproduce the correct relic density is typically smaller than in the former section. 
The exact value for A^r is now very dependent on the rest of the parameters. 

The input parameters for two explicit examples of this kind are given in Table 1 
and labelled aal) and aa2). We have chosen two cases with a pseudoscalar mass of 
m^o ^ 7 and 12 GeV, respectively. The values of Aat and Ax^^ for which the correct relic 
density can be obtained are represented in Fig. 6 for the range of RH sneutrino masses 
m^_^ = 6 — 45 GeV. We clearly see the above mentioned decrease in Aat when comparing 
these plots with the examples of the previous section. Regarding the CP-even Higgs 
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Figure 7: Diagrams contributing to the annihilation of RH sneutrinos into A^A^. 

spectrum, the lightest Higgs is SM-hke in both examples, with a predominant Hu 
composition and a mass of approximately 114 GeV. 

2.3 NN NN 

Let us finally address a last possibility that is characteristic of this model, namely 
annihilation into a pair of RH neutrinos, NN — )■ NN . Remember in this sense that the 
RH neutrino mass in this model is a consequence of radiative Electroweak symmetry- 
breaking and is thus related to the VEV of the singlet field as 

Mm = 2\mv, = 2^ . (2.15) 
A 

It is therefore possible to obtain a small value of the RH neutrino mass independently 
of the sneutrino mass (2.4), for which we still have two more free parameters to play 
with, namely the soft mass and the trilinear parameter Axj^- 

Since there is a lower bound in the value of the /i parameter in order to satisfy the 
experimental constraint on the chargino mass, > 105 GeV, and an upper constraint 
A < 0.6 in order to avoid Landau poles in the ROE of the NMSSM, a reduction in Mjv 
necessarily implies a small value of X^. For example, adopting the above constraints, 
one sees that Ajv < 0.018 (0.11) for M^v = 7 (40) GeV. 

The annihilation into a pair of RH neutrinos can proceed through s-channel Higgs 
exchange or t, w-channel neutralino exchange, as illustrated in Fig. 7. Notice that both 
the Higgs-neutrino-neutrino coupling and the sneutrino-neutrino-neutralino coupling 
are proportional to Aat {Cmmh° = —i/y/2XNS^o and C^^xi ~ —i\/2XNS^oN^_^). 
Therefore the upper bound on this quantity that we derived in the previous para- 
graph is a serious handicap in order to obtain a sufficiently large annihilation cross 
section. One possibility is to consider a resonant enhancement of the s-channel when 
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2m^^ rrifjo. However, if we want this channel to dominate over A^A^ — j- //, the 
corresponding Higgs exchange for the latter channel (which is obviously also enhanced 
by the resonance) has to be suppressed. The ratio of the contribution from both s- 
channels in the vicinity of the resonance with the lightest Higgs can be expressed as 
the fraction of the corresponding couplings, 




where the factor 6 takes into account the color factor of the hb channel and a factor 
1/2 for identical particles in the final state of the A^A^ diagram. Thus if we demand the 
above ratio to be larger than one, the lightest Higgs needs to be almost a pure singlino 
and its Hd composition extremely small. For example, for tan /3 = 3 and = 7 GeV 
it implies \SlolS]j,\ > 10 (30) for i?^^/,^ = 1 (10). 

This can be considered as a condition on the parameters entering the CP-even 
Higgs mass matrix (see e.g., expression (2.8) in [32]). In particular, the NMSSM input 
parameters can be chosen in such a way that the Ai's23 element of the Higgs mass 
matrix is larger than A^lis- In this sense, a small value of Ax is welcome but also a 
careful choice of the parameter k. 

The condition above clearly favours the use of small values of tan/3. This can be 
problematic since the supersymmetric corrections to BR(6 57) become sizable and 
it easily exceeds the experimental bound. Once more, avoiding this constraint further 
limits the choice of initial parameters. Similarly, the resulting value of a^^^ is very 
small and difficult to reconcile with experimental results from e+e~ data. Since the 
upper bound on A^r is relaxed for larger M^r, the amount of fine-tuning which is needed 
in order to obtain predominant annihilation into a pair of RH neutrinos decreases as 
Mjv increases. We have found numerical solutions for > 7 GeV with the correct 
relic density and satisfying the condition RNN/bl ^ 10- 

Finally, the t, w-channel can also be enhanced if the neutralino mass is decreased, 
however, this is generally not sufficient in order to recover the correct relic abundance 
for very light sneutrinos. 

Two explicit examples are given in Table 1 which satisfy the conditions described 
above, and are labelled as nnl) and nn2) in Table 1. For each of them, a scan is 
performed in the parameters A, Aat, and in order to account for all possible 
RH sneutrino masses and couplings. The value of the A and A^r parameters is set by 
relation (2.15), which we use to fix the RH neutrino masses to = 8 GeV in nnl) and 
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Figure 8: The same as in the right-hand side of Fig. 3 but for cases nnl) and nn2) of Table 1 
on the left and right, respectively. 

15 GeV in nn2). As it was already explained, the value of tan /3 is chosen to be as small 
as possible and the NMSSM parameters Ax, and k have also been fixed to values 
which lead to a small Ha component for the lightest Higgs. In particular, the lightest 
Higgs in case nnl) has a mass of approximately 50 — 60 GeV and \Sjjo\ ~ 0.004 — 0.01, 
|S'^o| ~ 0.2 — 0.4. Similarly, the Higgs in case nn2) has a mass in the range 70 — 90 GeV 
and \S^o\ ~ 0.01, l^^ol ~ —0.4 (due to the variation in A the mass of the Higgs also 
varies and these conditions can be fulfilled by a small range of masses). Since the 
correct relic abundance is only obtained when the resonant condition with the lightest 
Higgs is satisfied, this implies a value of the RH sneutrino masses of m^y^ ~ 25 — 30 GeV 
and 30 — 45 GeV, in cases nnl) and nn2), respectively. 

As explained above, if we consider lighter RH neutrino masses, the value of A^v 
has to be decreased. In our scans we have found that the correct RH sneutrino relic 
abundance could be obtained for Aat ~ 0.01 in case nnl) and A^v ~ 0.02 in case nn2). 
The viable points are represented in Fig. 8. Notice that these values of the coupling 
constant are considerably smaller than those obtained for the NN — > // case. 
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Figure 9: Diagram contributing to the spin-independent elastic scattering of RH sneutrino 
off quarks. 

3 Direct detection 

Let us now address the detectability of these particles in direct detection experiments. 
In general, WIMPs could be observed through their elastic scattering off nuclei (see 
Ref. [95] for a recent review), their interaction with quarks being described by an ef- 
fective Lagrangian that is valid in the non-relativistic regime where the collision takes 
place. In the case of RH sneutrinos there is only one Feynman diagram contributing 
at tree level to this process, namely, the t-channel exchange of neutral Higgses shown 
in Fig. 9. This leads to a Lagrangian describing the four-field interaction which only 
contains a scalar coupling, 

£e// 3 dq^NNm, (3.17) 

with 



3 



y ""T (3.18) 

where C^o^^ is the sneutrino-sneutrino-Higgs coupling, Yg. is the corresponding quark 
Yukawa coupling, and i labels up-type quarks {i = 1) and down-type quarks {i = 2). 
The spin-independent part of the sneutrino-nucleon elastic scattering cross section thus 
reads 

where rrip is the proton mass and 

— = E fk— + ^fTcJ2 —■ (3-20) 

qi=u,d,s qi=c,b,t ^* 

The hadronic matrix elements, fTq{= frq — fTq) and fxci— fro ~ Itg), are defined 
as {p\mgqq\p) = rripfi^^ and f^^ = 1 — X]g=tids /t?' determined experimentally as 
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/^^ = 0.020, /^^ = 0.026 and /^^ = 0.229. Being a scalar field, the effective Lagrangian 
contains no axial-vector coupling and thus the spin-dependent cross section vanishes. 
Using the explicit expressions of the quark Yukawa couplings, this quantity can be 
expressed as 

- 2MH.COS/3 ^ y"^^ V^' + + ~2r ) + ^ V"^" 27 

(3.21) 

where the term proportional to S^o corresponds to the interaction with the down-type 
quarks (the dominant contribution is due to the quark s) and the term proportional 
to S^o corresponds to up-type quarks. 

If the Higgs spectrum features a lightest Higgs with SM-like properties, i.e., with 
a mass of order 114 — 120 GeV, then the approximation 4m|^ -C m^o holds for very 
light sneutrinos and the above equation can be approximated as 

0.319 /p^ldli^y ,3,22) 



THp 2Mw COS f3 \ tan /3 

If, on the other hand, the lightest Higgs is lighter than the SM-like one (and neces- 
sarily featuring a larger singlet composition), the contribution from this lightest Higgs 
generally dominates and expression (3.21) can be approximated as 

A. ^ 0-31^ (v, + ^U,] ( . (3.23) 
nip 2Mw cos P \ tan /3 J y m^o J 

where we have used the quantities Vi and Ui defined in Eq.(2.10). Notice that the 
approximation 4m^ ^ m^o is not necessarily good now, since the Higgs can be very 
light, and this leads to the inclusion of the last factor. 

For moderate and large values of tan /3 the contribution from the s quark is the 
leading one to the spin-independent cross section. However, the second term can 
become sizable and even dominate for small tan /3 and if the lightest Higgs is mostly 
H^. Both in the case of the relic density, as well as in the scattering cross section, the 
condition that determines when the coupling from down-type quarks dominates over 
the coupling from up-type quarks is approximately the same, Vi ^ (0.4/ tan/3) Ui (see 
equations (2.13) and (2.14)). It should be emphasized that this rarely happens, but 
we include this possibility here for completeness. 

The resulting spin-independent contribution to the RH sneutrino elastic scattering 
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cross section off nuclei is then approximated as 



gj X f 0.31^ ^ , 0-42 ^ ^ 



"^^P TT (mp + m^J2 (^2MvK cos/3 J V ^ tan/3 ' 

In those cases where the hghtest Higgs is hghter than the SM one, a better approxi- 
mation (that incorporates the effect of the resonance) is 



^^P^ ^(^^ + ^^J2 (^2Mm/cos/3J + tan/3 \ ' ^^'^^^ 



We will use these expressions to extract some analytical predictions for the detectability 
of sneutrinos. Once more, in our numerical calculations the full expression for a^^^ has 
been included, without using any numerical approximations. 



3.1 NN ff 

Let us address first the case in which sneutrino annihilation into a pair of fermions is 
dominant. Using equations (2.8) and (2.10) the quantities V and U can be determined 
as a function of the lightest Higgs mass. Inserting these in Eq. (3.24) we obtain the 
following predictions for 

CT 0.13 rnittf^i 0.42 

2W^' if ^»r-^"' (3.26) 

ivp ^^(^^^ _|_ ^^j^^^z^ tan/3 



and 



,si 0.13 m^^m>,, P « 

m^(mp + m^J2(^,72^ — m2)3/2 ' tan/3 



In those cases where the lightest Higgs mass is small we can use expression (3.25) and 
substitute the quantities Vi and Ui to obtain 

m'l{mp + mfj^Y{m^^^— mlY/'^ \ m^o J ' tan/3 ' 

(3.28) 



and 



m^(mp + m^J2(^2^ _ ,77,2^3/2 I ,^2^^ / ' tan/3 

(3.29) 
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Figure 10: Spin-independent part of the RH sneutrino-quark elastic cross section as a 
function of the RH sneutrino mass for cases ffl) and ff2) and a scan in m^, Aat and 
Axj^. The black points correspond to those with the correct relic abundance, whereas the 
gray ones are those in which fi^/i^ is smaller than the WMAP constraint. The regions 
compatible with the CoGeNT observation of an irreducible excess and annual modulation 
are shown as the large (green) and small (red) areas bounded by a solid line. The shaded 
(blue) area is consistent with the DAMA/LIBRA experiment if no channelling effects are 
considered. The dashed, solid and dot-dashed lines correspond to the exclusion regions 
from the SIMPLE, CDMS and XENON experiments. 



In both cases there exists a correlation between the scattering cross section and the 
annihilation cross section but with a slightly different proportionality factor. Interest- 
ingly, imposing the correct relic density, and therefore the results for a^,^ and Ocg of 
equations (2.13) and (2.14), in the expressions above leads to a prediction of the cross 
section of order 10~^ pb for masses of order ~ 8 GeV, which can be compatible 
with the CoGeNT result. 

The theoretical predictions for are plotted as a function of the RH sneutrino 
mass in Fig. 10 for cases ffl) and ff2) of Table 1 for which we have performed a random 
scan in the m^y, Xj^ and Axj^ parameters, retaining only those points for which the RH 
sneutrino relic abundance is in agreement with the WMAP constraint (black dots) or 
smaller (gray dots). In case ffl) we have also represented the analytical approximation 
of Eq. (3.28) as a thick solid line, which turns out to be a qualitatively good approxi- 
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mation. Deviations happen because the numerical computation of the rehc density is 
more precise and also due to the contribution of the second-lightest Higgs to a^^. For 
illustrative purposes we also plot with a thick dashed line the theoretical predictions for 
the cross section if the term acc had been the dominant contribution to the annihilation 
cross section, as computed in Eq. (3.29). Similarly, in case ff2) the thick and solid lines 
correspond to the approximations Eqs. (3.26) and (3.27) obtained when annihilation 
into bb or into cc is dominant and including the contribution from the three CP-even 
Higgs bosons. The green area in both figures corresponds to the region consistent with 
the first CoGeNT results and the narrower red area is compatible with their latest 
ones. The cyan region is compatible with DAMA/LIBRA. Finally, the dashed, solid 
and dot-dashed lines correspond to the exclusion regions from the SIMPLE, CDMS 
and XENON experiments, respectively. 

Uncertainties in the hadronic matrix elements have not been included in the plot, 
but their effect is very easy to understand. The largest effect is due to the indetermi- 
nation in the strange quark content of the quark. This propagates into the theoretical 
predictions for and can be responsible for a variation of about an order of magni- 
tude [83, 84, 85, 86]. 

From our results we conclude that the predicted can be in agreement with 
the CoGeNT region without having demanded any further constraint, and solely as 
a consequence of the correlation between the diagrams that contribute to the RH 
sneutrino annihilation cross section and those for direct detection. This is however not 
only true for these particles, as it also happens with some other well-motivated WIMPs. 
In particular, this has already been pointed out for very light neutralinos both in the 
MSSM [21, 22, 25, 26] (although they are very constrained by low-energy observables 
[28]) as well as in the NMSSM [34]. A question worth investigating is then whether or 
not very light sneutrino dark matter can be distinguishable from other possible WIMPs. 
In this sense, including information from other sources of dark matter detection, such 
as indirect and collider searches, can shed some light on the specific model. 

3.2 NN 

Let us describe now the situation when annihilation into a pair of very light pseu- 
doscalars is dominant. Following from the discussion about the relic abundance in 
Section 2 we see how in this case there is no correlation between the diagram con- 
tributing to the direct detection (which is mostly dependent on the Xn parameter) 
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Figure 11: The same as in Fig. 10 but for cases aal) and aa2) in Table 1. 



and those for annihilation cross section (which are now dependent on a combination of 
various parameters). 

The imphcations for direct detection are easy to understand. If the annihilation 
cross section is dominated by the channel A^A^ — > A^A^, in the regions with the correct 
relic abundance the contribution from A^A^ — // is necessarily smaller than the one 
determined in the previous section. In other words, a^^i and Ocg in expressions (2.13) and 
(2.14) have now a lower value since the corresponding term a^o^o would provide the 
leading contribution to the calculation of the relic abundance. This entails a decrease 
in the predicted cx^^, according to Eqs. (3.28) and (3.29). 

As a consequence, points obtained in this scenario do not reproduce the CoGeNT 
results. Notice however that if the exclusion regions set by other experiments are taken 
at face value, these RH sneutrinos would still be allowed and survive as very light WIMP 
dark matter. We illustrate this in Fig. 11, where the theoretical predictions for as 
a function of the RH sneutrino mass are represented for cases aal) and aa2). As we 
can see, these very light WIMPs can have a scattering cross section which is not in 
tension with current exclusion regions. The predictions can be very low as in case aal) 
or within the reach of future experiments such as SuperCDMS. 

A similar phenomenon happens with very light neutralinos in the NMSSM when 
the resonant annihilation through a light pseudoscalar Higgs is invoked in order to 
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account for the correct relic abundance. In that case the predictions for sneutrino direct 
detection can be as low as ~ 10"^ — 10~^° pb [34, 38, 40]. Thus, although very light 
RH sneutrinos with annihilation into a pseudoscalar pair are clearly distinguishable 
from MSSM neutralinos, they might still be confused with neutralino dark matter in 
the NMSSM. However, we will see in the next section that the signals in colliders might 
differ. 



3.3 NN NN 

Let us finally address the scenario in which annihilation into a pair of RH neutrinos 
dominates. As in the case of annihilation into a pseudoscalar pair, the smallness of the 
Xn parameter in the regions with the correct relic density implies that the resulting 
spin-independent RH sneutrino-proton cross section is significantly suppressed. Once 
more, this happens because the Feynman diagrams that contribute to RH sneutrino 
annihilation are unrelated to those for direct detection. As a consequence, points of the 
parameter space where A^A^ — )■ A^A^ is the main annihilation channel would not account 
for the experimental results of the CoGeNT collaboration, if these are confirmed, but 
would still survive the bounds imposed by CDMS and XENON. 

This is clearly illustrated in Fig. 12, where the theoretical predictions for as 
a function of the sneutrino mass are represented for examples nnl) and nn2) of Ta- 
ble 1. The points with correct relic abundance accumulate in the narrow region where 
resonant annihilation is possible and can be as low as ~ 10"^ — 10^^*^ pb. 

We have also found it very difficult to obtain viable RH sneutrinos with m^j^ < 
20 GeV, as the parameters become extremely fine-tuned. From the point of view of 
direct detection, this example looks very similar to the one described in the previous 
section, and therefore also to very light neutralinos in the NMSSM. 



4 Invisible Higgs decay 

A common feature of most dark matter models featuring light WIMPs is the occurrence 
of a new invisible channel for the Higgs decay, namely the production of a dark matter 
pair. This also occurs in our scenario, since the light RH sneutrino has a sizable 
coupling to the CP-even Higgs bosons, and the resulting phenomenology can be very 
dependent on the specific scenario for sneutrino annihilation. 
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Figure 12: The same as in Fig. 10 but for cases nnl) and nn2) in Table 1. 

For any of the CP-even Higgses, the decay width of the process Hf — A^A^ can be 
expressed as 



T{H^ NN) 





2 

— m- 


4 





1 



for m^o > 2m^^ . (4.30) 



As in other models for very hght WIMPs, this decay mode can dominate for the 
hghtest Higgs. This can be the case of the very hght neutrahno, both in the MSSM 
and NMSSM, but also for a generic light singlet scalar [87, 88, 89, 90, 91, 92, 93, 94]. 
In our case the couplings C^oj^fj are very dependent on the properties of each of the 
annihilation scenarios that we presented in the previous Section, therefore we study 
each case separately. 



4.1 NN ff 

In order to estimate the invisible Higgs branching ratio we have to compare its decay 
width into RH sneutrinos with the decay width into fermions, which normally account 
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for tlie main visible decay cliannels'^. The ratio between the hghtest Higgs decay widths 
into RH sneutrinos and bb can be expressed as 

and we can define a similar quantity to compare with decays into cc, 

where we have used the approximation that m^„o ^ 4m^- . In Section. 2 we introduced 
two possible regimes in which the correct relic density could be obtained for a sneutrino 
annihilating into //. Depending on the lightest Higgs composition, the quantities P 
or U, defined in Eq.(2.9), provide the leading term for the annihilation into bb or cc, 
respectively. 

If we now impose that these sneutrinos reproduce the results from the CoGeNT 
experiment {a^ ~ 10~^pb ~ 2.6 x 10"^'^ GeV"^), the value of the coupling |C^ojv^| 
can be determined through equation (3.25), resulting in 

fl^^,,, « 2.8 X 10-' °^ ^ ' . (4.33) 

It can easily be seen that this ratio is typically much larger than one for realistic 
examples. For example, if the Higgs mass is small then the Higgs mass has to be 
mostly singlet, in which case \S^o\ becomes very small too and results in Rj^j^/^i ^ 1- 
If the Higgs mass increases, the composition becomes larger (closer to 1) and it 
can be explicitly checked that in the limiting case of a SM-like Higgs with a mass of 
order 114 GeV this ratio is still large. Remember in this sense that tan /3 is small in our 
scenarios. This already implies that the lightest Higgs, irrespectively of its composition 
and mass, has a very large invisible decay width. 

If we were in the regime where U ^ ^^T^ then Higgs decay into a cc pair would 
become larger than in bb (in fact, this coincides with condition under which sneutrino 
annihilation proceeds into bb or cc). In such a case, the ratio of this decay mode with 
Higgs decay into sneutrinos can be written as 

, 2 



sm 



(3 m-lo (mp + m^J' 



Rm/c, ^ 1-7 X 10-2 ' . , (4.34) 



•^If the decay into a pair of very light pseudoscalars, iJ" ^i^i, is open, its contribution can also 
be sizable. 
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Figure 13: Branching ratios of the decays of the lightest Higgs (left column) and second 
lightest Higgs (right column) in the possible different final states for cases ffl) and ff2) of 
Table 1. Only the points with correct relic abundance are plotted. 

and the same considerations as above would apply. 

Thus, the lightest Higgs will always tend to decay mostly into invisible particles. 
Notice that these results are independent on whether V or U dominates in Eq.(2.9). 

Regarding the decay of the second lightest Higgs, the RH sneutrino coupling to H2 
is less constrained by the relic abundance condition and therefore we can have several 
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possibilities. If the lightest Higgs is very light (e.g., below 100 GeV) and therefore 
singlet-like, the most conventional situation is that the second lightest Higgs is similar 
to that of the SM, with a mass of order 110 — 120 GeV, a large composition and 
a small but non- negligible Hd- Such a Higgs will decay mostly into a hb pair but it is 
likely that the presence of the very light sneutrino also induces a significant invisible 
decay width. 

This is actually the case of our example ffl). The branching ratios of the lightest and 
second-lightest Higgs are represented on the top row of Fig 13. Black dots correspond 
to the decays into a pair of RH sneutrinos and the rest of the decay products are 
indicated by different colours. As explained above is invisible, decaying mostly into 
RH sneutrinos. On the other hand, H2 (with a mass of 119 GeV) could be observed 
through its decay into hh. Notice however that this is significantly reduced with respect 
to an ordinary SM Higgs. 

A small change in the input parameters can alter significantly the phenomenology of 
the second-lightest Higgs. For example, in case ffl) if we take A\ = 500 GeV instead of 
550 GeV the invisible branching ratio increases considerably and can even be dominant 
for light RH sneutrino masses. This is shown in Fig. 14. 

Another possibility for the second-lightest Higgs is that, if the lightest Higgs is 
SM-like and with a mass around 110 — 120 GeV, H2 can be mostly singlet. In such 
a case, not only the light SM-like Higgs is invisible (as already explained above), but 
also this second Higgs. Case ff2) is one explicit example of this kind of scenarios. The 
branching ratios for and H2 for this example are depicted at the bottom row of 
Fig 13. 

Notice that in both cases the lightest pseudoscalar can also be relatively light, 
opening the decay channel — > A^A^. Although this channel can be very efficient, it 
is nevertheless not sufficient to compete with decays into RH sneutrinos. 

We therefore observe that, if the result of CoGeNT is imposed, the lightest Higgs 
is invisible, irrespectively of whether it is mostly singlino (and light) or the SM-like 
Higgs with a mass of around 110 to 120 GeV. Furthermore, in some scenarios also the 
second lightest Higgs can have a sizable invisible decay width. This is a very inter- 
esting property of this scenario and provides some potential discrimination criterium 
to distinguish it from the case of very light neutralinos. In particular, in Ref. [38] it 
was shown that in the case of very light NMSSM neutralinos the SM-like Higgs decay 
predominantly into a pair of very light Higgses. 
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Figure 14: The same as in Fig. 13 but with Ax = 500 GeV 

4.2 NN 

This example is qualitatively different from the one described in the previous Subsec- 
tion. In particular, the lightest Higgs is SM-like and therefore the annihilation into bb 
is typically large. The coupling of the lightest Higgs to the RH sneutrino pair is very 
dependent on the parameter space and this has an impact on the predictions for the 
invisible branching ratio. Now, since the coupling Xn takes smaller values than in the 
previous section we may expect that the contribution to the invisible decay width of 
all the CP-even Higgses is suppressed. 

This happens in case aal), where the lightest Higgs has a large component. 
Annihilation into bb constitutes now the main decay channel, thereby making it look 
similar to a SM Higgs. The predictions for the different branching ratios are displayed 
on the upper left-hand side of Fig. 15, where we observe that BR{Hi — A^A^) is 
significantly smaller than in the previous section. On the other hand, in case aa2) the 
first and second lightest Higgses are closer in mass and display a larger mixing. For 
example, despite having the same mass as in the previous example (114 GeV), the 
lightest Higgs has now a much larger singlet component. This is enough to enhance 
the predictions for the invisible branching ratio and, as we can observe in Fig. 15, 
BR{Hi —7- A^A^) > BR{Hi — )■ bb) for the whole range of RH sneutrino masses. 
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Figure 15: The same as in Fig. 13 but for cases aal) and aa2) of Table 1. 



The second-lightest Higgs is normally singlet-like and as a consequence its anni- 
hilation into SM particles is already suppressed with respect to annihilation into RH 
sneutrinos and we should expect a sizable BR{H^ — > A^A^) as well as BR{H^ — > 
from the annihilation into a pair of very light (singlet-like) pseudoscalars. This is in- 
deed what happens, as displayed on the right-hand side of Fig. 16 for cases aal) and 
aa2) on the top and bottom, respectively. Again, the exact values are very dependent 
on the specific choice of the input parameters. 
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Figure 16: The same as in Fig. 13 but for case nnl) of Table 1. 

This situation differs from the case of very hght neutrahnos in the NMSSM. In 
the regions where the neutrahno has the correct rehc density due to the annihilation 
into a pair of very hght pseudoscalars, the SM-hke Higgs annihilates preferentially into 
a pair of pseudoscalars Ref. [38]. However, in our case we observe that this is not 
necessarily the case and the SM-like Higgs (which in our examples coincides with H^) 
either appears as a Higgs with SM-like decays (into bb) or with a significant invisible 
decay. 

4.3 NN NN 

The Higgs spectrum in these examples features a light scalar Higgs, with a large singlet 
composition, and a second-lightest Higgs which is SM-like. In fact, as we stressed in 
Section 2.3, the composition of the lightest Higgs has been carefully chosen in such a 
way that it couples preferentially to a pair of RH neutrinos (so that annihilation into 
these is the main channel). More specifically, we chose the input parameters in such a 
way that 5*^0, 5'^o ^ S^o is fulfilled. It is then clear that the main decay modes of the 
lightest Higgs are — > A^A^ (if 2m^^ < mj^o) and — ?■ A^A^. Notice in this sense 
that while the decay into RH sneutrinos is an invisible channel, this needs not be the 
case for the decay into RH neutrinos. 
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In fact, RH neutrinos in our model are not stable particles. They couple to a 
lepton and a boson through their mixing with LH neutrinos, and thus can undergo 
a three body decay IW* -> IIul, N IW* Iqq, N vlZ* ulU or 

— )■ ulZ* — )■ uigq. Since the RH neutrino coupling to a lepton and a W boson 
is suppressed by the small Yukawa coupling (which determines the neutrino LH-RH 
mixing) the lifetime of RH neutrinos can be either long enough for them to produce 
displaced vertices or too large, thereby counting also as an invisible component. 

If Fig. 16 we represent the resulting branching ratios for the different decay modes 
of the lightest and second-lightest Higgses in example nnl). Black dots correspond 
to BK{Hi — )■ NN), whereas empty circles represent BR(ifJ' — NN). As commented 
above, these constitute the main decay modes for the lightest Higgs. 

Since the value of the coupling Aat is very small in these examples, the coupling of 
the RH neutrino and sneutrino to the heavier Higgses (which are, respectively, mostly 
Hu and Hd) are small. Thus there is no contribution from these modes to the branching 
ratios of these. Still, since we are dealing with the NMSSM, the SM-like Higgs can 
have exotic decays, although these are more dependent on the particular choice of input 
parameters. For example, in case nnl) the channel H2 — ?■ H^H^ is open and provides 
the main contribution to the decay width of Also, the neutralino is a relatively 
light singlino (with a mass of 40 GeV) and also provides a prominent decay mode. 
In sum, in this last case the second lightest Higgs (which is SM-like) can have a very 
similar phenomenology to the case of neutralino dark matter in the NMSSM. 

Summarising, this scenario could potentially be discriminated from other models 
with very light WIMPs through the study of the displaced vertices produced after the 
production of a pair of RH neutrinos and their subsequent decays. The signature of 
this decay mode at the LHC depends largely on properties of the RH neutrino which 
are intimately related to the details of the see-saw mechanism. This information is 
contained in the structure of the neutrino Yukawa coupling. In particular, the RH 
neutrino lifetime, its decay products (especially whether it decays "democratically" in 
electrons, muons or taus) and the effect of having three RH neutrinos (and sneutrinos) 
have to be carefully investigated. This analysis is beyond the scope of the present work 
and will be the subject of a more specific future study. 
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5 Indirect detection 



As in most of the WIMP scenarios, annihilations of very hght RH sneutrinos in the 
Galactic halo may contribute to the observed fluxes of high energy p/p, e+/e~, v/i? 
and gamma rays. Among them, we will here address the possible signatures of our 
model in the gamma ray flux from the Galactic Centre, a region which is currently 
being observed by the Large Area Telescope on the Fermi satellite (Fermi-LAT) [96]. ^ 

In the usual neutralino dark matter scenario, annihilation in the Galactic halo tends 
to be suppressed due to the p-wave dominance of the annihilation cross section because 
of the Majorana nature of the neutralino (see, e.g., the recent study of the gamma 
ray flux for neutralinos in the cMSSM in Ref.[99]). On the contrary, RH sneutrino 
annihilation is s-wave dominant, thus we expect a larger gamma ray flux. 

As explained in Section 2, the dominant channel for sneutrino annihilation can be 
//, A\A\ or A^A^, depending on the chosen parameters. Among them, A^A^ mode dom- 
inant scenarios are potentially interesting because the gamma ray spectrum produced 
in the subsequent decay of the RH neutrino into three fermions through Electroweak 
interaction (A^ — Uvli N — > Iqq oi N ^ QQ^l) is distinctive from the annihilation 
products of the // channel. This possibility is not present in other WIMP models and 
thus is not generally addressed in the literature. 

In general, the photon spectrum from WIMP annihilation or decay has three com- 
ponents, the prompt, the inverse Compton scattering (ICS) and synchrotron radiation. 
The prompt component involves both continuous and line emission. The main sources 
of continuous spectra are 27 decays of 7r° and final state radiation. To calculate the 
expected spectrum we have written a routine generating events of three body decays 
of A^ with an appropriate weight using the PY3ENT routine in PYTHIA 6.4 [100]. For 
the hb and A^A^ modes, we also used PYTHIA 6.4. In Figs. 17 and 18, we present the 
energy spectra of gamma rays produced in the subsequent decays of final state particles 
for each annihilation mode. Although in various WIMP models the annihilation into 
27 appears at the one-loop level, it is considered as a smoking gun channel of indirect 
dark matter detection because it is distinguishable from any conceivable astrophysical 
background. However, for very light RH sneutrino dark matter, the line emission tends 
to be buried in the large astrophysical background after taking into account the energy 
resolution of the apparatus. Thus we can neglect its contribution in the analysis. 

^Although we do not discuss here the possible signals in the high energy cosmic ray (CR) like e+e~ 
or pp, these may also constrain light dark matter candidates [97, 98] . 
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Figure 17: Left: dN^/dE from NN bb. Right: dN^/dE from NN 
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Figure 18: Left: dN^/dE from NN NN. Right: dN^/dE from A^A^ ^ A\A\ 



Other contributions are produced from the propagation of charged stable prod- 
ucts hke e+/e~ and p/p. Once they are produced, they are diffused and loose their 
energies due to ICS with ambient start light or cosmic microwave background ra- 
diation, and synchrotron radiation by the Galactic magnetic field. Then the same 
energy loss mechanism of charged cosmic ray also generates photons. For heavy 



37 



enough dark matter, these contributions can exceed the prompt one in the relatively 
lower energy region [101]. However, for light dark matter with mass ~ (9(10GeV) 
the maximum energy of gamma rays from ICS by starlight with energy u ~ leV 
is around [E/meYu ~ O.lGeV, which is near the threshold energy of Fermi-LAT. 
Synchrotron radiation on the other hand is produced in the radio frequency range 
around v,y,^ - (1/3) 3eBE^/47rml = 1.3 kHz (B/lmG) {Ejm^f < 10^ GHz (5/lmG) 
[101, 102]. Since the constraints arising from the observation of radio frequency are be- 
yond the scope of this article, in the present analysis we neglect the contributions com- 
ing from both ICS and synchrotron. Being able to ignore the ICS and the synchrotron 
radiation makes the analysis much simpler because of the absence of complications 
stemming from uncertainties in CR propagation model. 

Following the standard procedure [82] the gamma ray flux can be written as 



(Afi) = > / dV 

dE"- ' ^Jn<An 2 Ans^ dE 



El 



r n'^(r) dm 

/ ds dn / ^ ' J aiv)—^ 



^^j:ia..)^J(Amn, (5.35) 



''TV, 



where Ubm is the dark matter number density, i represents all the possible annihilation 
channels {ee, fifi, rf, uu, dd, cc, bb, A^A^, NN} , dN^/dE{E) is the expected number 
of photons in the energy range of {E, E + dE) produced from a given process i, and 
the so-called halo factor J(Af2) is defined as 

J{fl)An = I dn I ds p2 ( r = Jrl + s"^ + 2rQS cos 9] . (5.36) 

Pe^Q J An Jlinc-of-sight V / 

Here we have used the canonical value Pq = p{rQ) = 0.3 GeV for the dark matter 
density around the Sun^ and a distance to the GC of = 8.5 kpc. The halo properties 
are thus factorized and encoded in the single factor J{AQ). In order to take into 
account the possible astrophysical uncertainties, we use two halo models which are 
supported by N-body simulations, namely NFW [108] and Einasto [109] as well as the 



^More recent determinations in Rcf. [103, 104, 105, 106, 107] indicate slightly larger values. In any 
case, the fluxes corresponding to different values of the local density can be easily deduced from the 
appropriate scaling. 
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9.6 X 10-^ 


2.4 X 10^2 


4.6 X 10-2 


9.6 X 10-2 


NFW 


1.35 


5.95 


7.91 


10.5 


Einasto 


2.10 


11.6 


15.2 


19.7 


Isothermal 


0.0130 


0.319 


0.615 


1.22 



Table 2: Halo factors multiplied by given solid angles, J(Af2)Ar2, for the different halo 
models, NFW, Einasto and isothermal model. ROI around GC with radius 1°, 5°, 7° and 
10° are used. 

isothermal halo model [110, 111] as a reference. The corresponding density profiles are 
parametrized as 

r PsTs/ril + r/rsf , NFW 
Pir) = l exp [-|(r" - 1)] , Einasto (5.37) 
y Ps/ {1 + {r /tsY) , Isothermal 

where = 20 kpc and 5 kpc for NFW and isothermal model, respectively, and a = 
0.17. The value of ps is fixed to reproduce the dark matter density around the solar 
system. In Table 2, we calculate J{AQ)AQ for each halo model considering a region 
of interest (ROI) around the GC with radius 1°, 5°, 7° and 10°. 

To compare the predicted gamma ray flux from RH sneutrino annihilation with 
the gamma ray flux currently observed by Fermi-LAT, we have used data from the 
Fermi Science Support Center (FSCC) archive [112] ^. We have selected the ROI as a 
circular region with radius 1°, 5°, 7° and 10° around the GC with RA= 266.46° and 
Dec= —28.97°. Following the suggestion of the FSCC, we used a zenith angle cut 105° 
to reject photons coming from the Earth. We used gtbin tool to make 20 bins which 
is equally spaced in logarithmic scale in energy from 0.1 GeV to 100 GeV. 

To obtain the flux map for gamma rays from the counts map of photons actu- 
ally detected in the experiment, we divided the count map by the exposure map in 
each position and energy bins. This can be calculated from the spacecraft data and 
the instrument response function using the gtltcube, gtexpcube2 tools. We used 

^We extracted gamma ray data taken from 4th of August in 2008 (15:43:37) to 4th of July in 
2011 (14:08:11). We selected signals classified as DIFFUSE only, which are appropriate to analyse the 
diffuse gamma ray emission using the gtselect tool in Fermi Science Tools [113]. In both processes 
of selecting data and calculating the exposure map, we used the specific instrument response function 
(IRF) P6_V11_DIFFUSE to be consistent. 
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the gtexpcube2 tool to reflect the azimuthal dependence of the effective area of the 
apparatus. 

As described in [114], we have taken the systematic uncertainty in the effective area 
of the apparatus as 10 % at 0.1 GeV, decreasing to 5 % at 0.56 GeV and increasing to 
20 % at 10 GeV, and interpolated between the points. We have extrapolated above 10 
GeV with a constant value. Since the systematic uncertainty in the background model 
is not fully understood yet, we used a fixed background emission template provided by 
the Fermi Science Team [115]. 

There are three main components of background^, the diffuse galactic emission 
(DGE), the resolved point sources (PS) and the isotropic gamma ray background (IGB). 
Here we used the DGE model map gll_iem_v02 . f it and IGB model supplied by the 
Fermi Science Team. For PS, we used the preliminary second Fermi-LAT catalog found 
in [116]. PS in the catalog are fitted in PowerLaw, LogParabola and PLExpCutoff 
whose form can be found in the accompanying draft. We modeled its contribution by 
summing over the modeled fluxes of PS in a given ROI. 
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m^v (GeV) 


rrij^o (GeV) 
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3.47 X 10^26 


bb 
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3.47 X 10^26 


cc 
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aal-10) 


4.52 X 10-26 




10 
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aal-20) 


3.19 X 10-26 
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aal-40) 


2.90 X 10-26 




40 
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aa2-20) 


3.19 X 10-26 
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11.7 


aa2-40) 


2.90 X 10-26 
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11.7 


nn-23) 


3 X 10-26 
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23 
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nn-30) 


3 X 10-26 
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30 
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3 X 10-26 
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40 


15 


* 


nnB-22) 


3 X 10-28 
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22 
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nnB-25) 


3 X 10-25 


NN 


25 


8 


* 



Table 3: Benchmark points chosen to calculate expected gamma ray fluxes. We assume 
sneutrino annihilates only into specified mode. The parameters marked with * are not 
relevant in the analysis. 

''Here we use the terminology background although we include both background and foreground in 
the analysis. 
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To illustrate the various annihilation channels and mass spectra in our model we 
have chosen 12 representative benchmark points as in Table 3, inspired in the examples 
discussed in the previous sections. For simplicity, we assumed that annihilations occur 
only through the given specific channels. 

5.1 NN ff 

Let us first consider the case in which RH sneutrinos annihilate into a pair of fermions. 
As we explained in Section 2, the dominant channel is either bb or cc, depending on the 
choice of parameters, while the rr channel is always negligible. We have chosen = 
8 GeV which is compatible with the CoGeNT result. The benchmarks where either 
annihilation into bb or cc dominates are thus labelled bb-8) and cc-8), respectively. 

The predicted gamma ray fiux is represented as a function of the energy in Fig. 
19 for the various choices of ROI. As we can observe, the fiux with Einasto or NFW 
profile is larger than the fiux from resolved point sources. However, the expected fiux 
calculated for 5°, 7° and 10° is still smaller than DGE, and thus still consistent with 
the observed data after taking the systematic uncertainties into account. 

Notice that the predicted fiux for 1° is comparable to the DGE component. How- 
ever, in this region, a significant discrepancy between fiuxes from background model 
and observed fiux is present. Indeed, our knowledge of the background in this region is 
still poor. In particular, the suppressed fiux observed at low energy cannot be explained 
by point sources near the GC which are not resolved yet. 

5.2 NN 

In the second set of benchmark points, we focus on examples in which the RH sneutrino 
annihilates into A^A^. We consider two cases aal) and aa2) defined in Section 2 and 
fix the RH sneutrino mass to m^y^ = 10, 20 and 40 GeV for aal), and =20 and 
40 GeV for aa2). The various benchmark points are labelled accordingly in Table. 3. 

The main difference in the predicted gamma ray spectra for scenarios aal) and aa2) 
is the mass of the lightest pseudoscalar Higgs. In cases aa2) the pseudoscalar is heavy 
enough to decay into bb, whereas in cases aal) this is not possible As shown in Fig. 
18, the presence of the bb mode leads to a softer spectrum. This is further evidenced 
in Figs. 20 and 21, where the predicted gamma ray fluxes for the different benchmark 
points are represented. Contrary to the // scenarios, we have varied the mass of RH 
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Figure 19: Expected gamma ray flux in bb-8) and cc-8). ROI with radii 1° (top left), 5° 
(top right), 7° (bottom left) and 10° (bottom right) are used. 

sneutrinos in the relatively large range from 10 GeV to 50 GeV in these cases. The 
number density of dark matter with a fixed halo model is inversely proportional to the 
RH sneutrino mass-squared, thus so is the gamma ray flux. 

These results are similar to those obtained in the // scenarios analysed previously. 
The expected gamma ray fluxes are too small to be observed except for the case where 
a ROI of 1° is considered. Therefore, without a significant improvement of our under- 
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Figure 20: Expected gamma ray flux in aal-10), aal-20) and aal-40). ROI with radii 1° 
(top left), 5° (top right), 7° (bottom left) and 10° (bottom right) are used. 

standing of the background, we cannot constrain the relevant parameter space of these 
examples with gamma ray flux from the GC. 
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Figure 21: Expected gamma ray flux in aa2-10), aa2-20) and aa2-40). ROI with radii 1° 
(top left), 5° (top right), T (bottom left) and 10° (bottom right) are used. 

5.3 iViV ^ iViV 



The last class of benchmark points corresponds to scenarios in which RH sneutrino into 
A^A^ is dominant. These cases are potentially very interesting since the RH neutrino 
subsequently decays into three fermions (A^ — )■ llvi^ or N ^ ^q)- This leads to a 
distinctive gamma ray spectrum. 
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Figure 22: The distributions of (de-)enhancement factor {(tv)x=o/ {o'v)x=x; in RH sneutrino 
mass for nnl) case (left panel) and nn2) case (right panel). 



Moreover, as explained in Section 2.1, the correct RH sneutrino relic abundance in 
these cases is only obtained through the CP-even Higgs resonance. In such a case, 
there can be a significant enhancement (or suppression) of its annihilation in the dark 
matter halo and consequently in the gamma ray flux. This so-called Breit-Wigner 
enhancement has been studied in various models and contexts [117, 118, 39]. 

To illustrate this effect we have represented in Fig. 22 the enhancement factors of 
the annihilation cross section, defined as {av)x=o/ {o'v)x=Xf for the points with the cor- 
rect RH sneutrino relic density in examples nnl) and nn2) of Table 1. As we observe, 
the enhancement factor ranges from 10~^ to 10^, depending on the difference between 
2m ff and mn, and the decay width of corresponding CP-even Higgs. Consequently, 
these benchmark points cannot be represented with a single fixed cross section. We 
have therefore defined cases nn-23), nn-30) and nn-40) with a canonical annihilation 
cross section 3 x 10~^^ cm^/s, while 3 x 10~^^ cm^/s and 3 x 10~^^ cm^/s are used 
for the nnB-22) and nnB-25) cases, respectively, to study the impact of the enhance- 
ment/suppression on the gamma ray flux. 

The corresponding predictions for the gamma ray flux are represented in Fig. 23 
for benchmark points nn-23), nn-30) and nn-40). We can see again that the expected 
gamma ray flux is too small to constrain RH sneutrino annihilation in the halo in these 
cases except for a ROI of 1°. However, these gamma ray fluxes are calculated under the 
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Figure 23: Expected gamma ray flux in nn-23), nn-30) and nn-40). ROI with radii 1° (top 
left), 5° (top right), 7° (bottom left) and 10° (bottom right) are used. 

assumption that RH sneutrinos annihilate with fixed annihilation cross section with 
value of 3 X 10~^^ cm^/s, while Fig. 22 shows that the annihilation cross section in the 
halo can be at least one order of magnitude larger or smaller. Therefore, the fluxes 
shown in Fig. 23 should be adjusted by the appropriate factors. 

The impact of this Breit-Wigner enhancement is illustrated in Fig. 24, where two 
cases with the suppression factor of 10~^ for nnB-22) and the enhancement factor of 10 
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Figure 24: Expected gamma ray flux in nn-23), nnB-22) and nnB-25). ROI with radii 1° 
(top left), 5° (top right), 7° (bottom left) and 10° (bottom right) are used. 

for nnB-25) are represented. With the Einasto model, the expected fluxes of nnB-25) 
are larger than the observed fluxes. In other words, the gamma ray flux from the GC 
region observed by Fermi-LAT is already constraining a portion of the parameter space 
of our model. Nevertheless, a thorough analysis requires scanning over full parameter 
space of the model, and we leave it for a future study. 
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6 Conclusions 



In this paper we have investigated the viabihty of very hght RH sneutrinos in the 
NMSSM and analysed the imphcations for direct dark matter detection, the poten- 
tial effects on Higgs phenomenology and the prospects for indirect detection through 
gamma rays. The model contains a new singlet superfield that provides RH neutrinos 
and sneutrinos and three input parameters, a soft scalar mass m^, a coupling constant 
Xn and the associated trilinear parameter Axj^. 

First we have studied in detail the conditions under which RH sneutrinos in the 
NMSSM can be very light and reproduce the correct value for the relic abundance. We 
have found that this is possible in three different scenarios, namely when RH sneutrino 
annihilates mainly in fermions {NN — )■ //), in a pair of very light pseudoscalar Higgses 
{NN A? A?), or in RH neutrinos {NN NN). 

The first case (A^A^ //) is possible in the presence of a light (singlet-like) scalar 
Higgs or with an increase of the coupling constant Aat. In both cases, the diagrams 
contributing to the sneutrino annihilation in the early Universe and those describing 
their scattering cross section off quarks are correlated. We observe that this correlation 
implies that if the correct relic density for RH sneutrinos is imposed, their predictions 
for direct detection are compatible with the WIMP interpretation of the CoGeNT 
results. We derive an analytical approximation to make this correlation explicit. This 
scenario has a deep impact in the predicted invisible decay width of the Higgs bosons. 
We show that the lightest Higgs is mostly invisible since it decays in a pair of RH 
sneutrinos. Interestingly, this applies to both a very light singlet-like Higgs and a 
SM-like Higgs with mass of order 114 GeV. In fact, we also observe that the second 
lightest Higgs in this scenario can also have a significant invisible decay width, although 
whether or not the invisible modes dominate is much more dependent on the specific 
input parameters. 

The second case (A^A^ — A^A^) is possible when the pseudoscalar Higgs is light 
and singlet-like. There is now no correlation between the annihilation cross section 
and scattering cross section off nuclei and as a consequence the theoretical predictions 
for RH sneutrino direct detection cannot account for the CoGeNT signal and are much 
smaller, of order 10^^ — 10^^° pb. However, they are consistent with the exclusion re- 
gions set by the CDMS and XENON experiments. In some cases they could be within 
he reach of future experiments such as Super CDMS. The resulting Higgs phenomenol- 
ogy is slightly different than in the previous scenario. The invisible decay width of the 
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lightest Higgs is not necessarily large for the lightest Higgs (which is SM-like), but it is 
typically sizable for the second lightest Higgs (although much more dependent on the 
specific input parameters of the model). 

The third possibility (A^A^ — A^A^) is more constrained but potentially very in- 
teresting. We described the stringent conditions under which this diagram can be 
dominant, which is only possible in the resonance with a light singlet-like Higgs with a 
very small Hd component. Once more there is no correlation between the annihilation 
cross section and scattering cross section off nuclei and the resulting a^^^ is very small, 
of order 10~^ — 10~^° for RH sneutrinos with masses ^20 GeV and featuring RH 
neutrinos with masses above Mat > 8 GeV. Although it cannot explain the CoGeNT re- 
sult, it would be compatible with the bounds set by CDMS and XENON. This scenario 
leaves a potentially characteristic signal in colliders, namely the lightest Higgs decay 
into a RH neutrino pair whose subsequent (late) decay leaves two displaced vertices. 

Last we have investigated the theoretical predictions for the gamma ray flux from 
the GC and compared it with the results from the Fermi satellite. In this respect, 
the cases A^A^ — )■ // and A^A^ — )■ A\A\ lead to the conventional results obtained for 
other DM models. Much more interesting is the new annihilation mode (A^A^ — j- NN), 
due to the subsequent decay of the RH neutrino. We have characterised this new 
possibility finding that the resulting spectrum is an intermediate case between those 
of annihilation into bb and rf . Furthermore, since the annihilation occurs at the Higgs 
resonance, a Breit-Wigner enhancement of the resulting annihilation cross section in the 
dark matter halo can result in a boost factor as large as a factor 100 (or a suppression 
factor of several orders of magnitude). 
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